Ionization of clusters in intense laser pulses through collective electron dynamics 
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The motion of electrons and ions in medium-sized rare gas clusters (~ 1000 atoms) exposed to intense 
laser pulses is studied microscopically by means of classical molecular dynamics using a hierarchical 
tree code. Pulse parameters for optimum ionization are found to be wavelength dependent. This 
resonant behavior is traced back to a collective electron oscillation inside the charged cluster. It 
is shown that this dynamics can be well described by a driven and damped harmonic oscillator 
allowing for a clear discrimination against other energy absorption mechanisms. 

PACS numbers: 36.40.Gk, 31.15.Qg, 36.40.Wa 



The interaction of of intense laser radiation with clus- 
ters has been of continuing interest Q, 0, IE 0, IE IE i 
pushed by prominent findings as the emission of keV- 
hotons Q , highly charged ions 0, Q , or fast fragments 
Basically, all these phenomena are caused by the 
exceedingly effective absorption of energy from the laser 
field into the cluster. This enhanced absorption (when 
compared to atoms or bulk matter) is possible due the 
initially solid-like atomic density in the cluster in combi- 
nation with the rapid expansion of the cluster on a fem- 
tosecond time scale, i. e. typically during the laser pulse. 
More detailed insight into the mechanism of energy ab- 
sorption can be gained by pulse length variation. Using 
pulses of equal energy, i. e. longer pulses have lower in- 
tensities, it has been found that there is an optimum 
pulse length with maximum absorption [E EJ • This holds 
true for small rare gas and metal clusters (of the order 
of some 10 atoms) as well as for large clusters, although 
the underlying reason for maximum absorption can be 
very different and is in fact an issue of current debate. 
For small rare gas clusters under laser pulses of peak in- 
tensities in the tunneling regime (/ > 10 15 W/cm 2 ) such 
an optimum could be attributed to the mechanism of en- 
hanced ionization Q known from molecules For small 
metal clusters exposed to similar laser pulses the exis- 
tence of optimal absorption was interpreted as a plasmon 
resonance phenomenon 6] in analogy to the well known 
dipole resonance of the valence electrons in perturbative 
photo-absorption ^(j or low-intensity laser pulses [Tl| . 
For large clusters (of more than 10 5 atoms) it has been 
proposed that strong laser pulses create a nanoplasma 
inside the cluster 0, In the course of the expan- 
sion of the cluster the electron density and consequently 
the plasma frequency decreases resulting in strong energy 
absorption at resonance with the laser. 

Separating the different mechanisms from each other 
requires specific and clear signatures for each process 
which are difficult to identify in a multi particle system 
such as a cluster. Surprisingly, dipole resonant absorp- 
tion dynamics in a cluster can be very well characterized 
by a simple driven damped harmonic oscillator, which 
describes the dipole response of the electrons inside the 
cluster. 

In the following we will demonstrate the validity of this 



simple description with full dynamical microscopic calcu- 
lations for Xenon clusters (~ 10 2 . . . 10 3 atoms) in strong 
optical laser pulses (wavelengths A = 520 ... 1170 nm, 
intensities / ~ 10 14 . . . 10 16 W/cm 2 , pulse lengths T ~ 
10 . . . 1000 fs). Our approach is similar to those used be- 
fore for intense laser-cluster interaction |l4|. However, 
we have been forced to use a completely new propagation 
scheme, namely a hierarchical tree code [l5j. to handle 
of the order of 10 4 charged particles (~ 1000 ions and 
^8000 electrons) with their mutual interactions. Orig- 
inally developed for gravitational TV-body problems in 
cosmology [Tg, the hierarchical tree code allows us to 
follow the dynamics of all charged particles over a few 
hundred femtoseconds with typical time steps of attosec- 
onds. 

Atoms are initially arranged in so-called Mackay icosa- 
hedra |17| known to be optimal structures of Lennard- 
Jones clusters. We distinguish between free and bound 
electrons, whereby the latter ones are not treated explic- 
itly. The condition for creation of an electron is that none 
of the other electrons has a negative binding energy to 
the ion under consideration. In this case a new electron 
is "born" at the position of that ion with a kinetic energy 
to satisfy the ionization potential. The charge of the ion 
is increased by one. Free electrons as well as ions are 
described classically as charged particles in an oscillating 
field interacting via a softened Coulomb interaction W. 
This is defined for a pair of particles with charges q\ and 
q-i and a distance of ryi as W = q\ 52/VV12 2 + 1. 

First of all, we present the calculated pulse length de- 
pendence for three different laser wavelengths A. Fig- 
ure n shows the final averaged charge state per atom 
from Xe56i clusters after laser impact as a function of 
the pulse length T or the peak intensity /, respectively. 
The applied laser pulse was linearly polarized with a field 



F(t) = J ^ I ^ COs2 (5t) COS ^ for 1*1 < T (i) 
elsewhere 

where la — 3.51T0 16 W/cm 2 . In order to keep the en- 
ergy of the different pulses constant we fixed the prod- 
uct of intensity and pulse length / x T = Iq x 4fs. For 
short pulses (T < 30 fs) , where the cluster atoms have not 
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FIG. 1: Average charge per atom from Xes6i clusters after 
laser impact according to Eq. Q as a function of the pulse 
length T or the peak intensity I, respectively, for three dif- 
ferent laser wavelengths A. The energy of the pulse, i. e. the 
product / x T, is kept constant. 



enough time to react on the charging, the final charge 
state decreases with an increasing pulse length due to the 
lower intensity. For longer pulses (T>50fs), however, 
the final charge state increases despite the intensities be- 
come smaller. This can only be understood if one consid- 
ers the expansion of the cluster, for a detailed explanation 
see below. Finally, for very long pulses (T > 400 fs) the 
cluster is already completely fragmented before the laser 
pulse reaches its peak intensity rendering the ionization 
similar compared to the case of single atoms. Qualita- 
tively, this behavior is the same for all three frequencies 
(Fig. nj. However, the shift of the optimal pulse length 
towards longer pulses for longer wavelengths is charac- 
teristic for a resonant ionization mechanism. 

To gain insight into the mechanism of ionization we 
discuss the dynamics of a Xeg23 cluster in a laser pulse 
with 7 = 9- 10 14 W/cm 2 , cf. Fig. [2 The pulse has short 
rise and fall times of 20 fs and a long plateau of 160 fs in 
order to eliminate effects from the time dependence of the 
laser pulse itself. As can be clearly seen in Fig. the 
cluster ionization occurs in two steps and the total cluster 
charge Q (sum of total ionic charge Qion and charge of 
those electrons which are inside the cluster volume) as a 
function of time t is well represented by 



Q(t) = E 



i=l,2 



l + cxp(-(t-U)/SU) 



(2) 



In the 1st step during the rising of the pulse (ti ± 5t\ = 
— 86fs±3fs) electrons are emitted mainly due to field 
ionization. This process slows down, however, already 
before the plateau intensity is reached due to the increas- 
ing space charge (cluster charge Q w 1000 at t = — 80 fs). 



Because of this attractive space charge, one may distin- 
guish between inner ionization, which accounts for excita- 
tion from localized electrons to quasi-free electrons mov- 
ing inside the cluster volume, and outer ionization, which 
corresponds to the final escape of the quasi-free electrons 
from the cluster into the continuum. Note, that the space 
charge can hold an appreciable number of quasi-free elec- 
trons which engage into collective motion discussed be- 
low, in contrast to the ionization dynamics of molecules 
or small clusters with almost no quasi- free electrons @. 
The charging up of the cluster leads to an expansion as 
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FIG. 2: Dynamics of Xeg23 in a strong laser pulse (A = 
780 nm, 7 = 9- 10 14 W/cm 2 , rise and fall time 20 fs, plateau 
for t — —80 ... + 80 fs). All quantities are shown as a func- 
tion of time t. a: Average charge per atom (circles and fit 
from Eq. j^). left axis) and corresponding rate (gray filled 
line, right axis), b: Radii R of all cluster shells in units of 
their initial radii Rq. c: Centre-of-mass velocity « C m of the 
electronic cloud inside the cluster volume. Note, that the 
oscillations are spatially along the linear polarization of the 
laser, whereas the electron velocity perpendicular to the laser 
polarization is very small and hardly to see in the figure, d: 
Phase shift <f>t of the collective oscillation in laser direction 
with respect to the driving laser, see text. 



can be seen for t > — 60 fs in Fig. [21). During the expan- 
sion a 2nd ionization step occurs which lasts for a much 
longer time and leads to an increase of the average ionic 
charge from about 1 to more than 6. During this time, 
the quasi-free electrons in the cluster are driven collec- 
tively back and forth along the polarization direction of 
the laser which is evident from their centre-of-mass (cm) 
velocity v CM shown in Fig. [2^. 

This oscillation can be modelled by a driven and 
damped classical harmonic oscillator 



X(t) + 2T t X(t) 



VL t 2 X{t) 



F(t) 



(3) 



with X(t) the CM position of the electron cloud, F(t) 
the driving laser amplitude, and £lt and Tt the eigen- 
frequency and damping rate, respectively, which are de- 
termined by the cluster. The index t indicates that 
due to ionization and expansion of the cluster, both, 
£lt and T t , may depend parametrically on time. Under 
the assumption of a spherical, uniformly charged clus- 
ter with total ionic charge Q lon and radius R the po- 
tential inside the cluster is harmonic with an eigenfre- 
quency £lt = \/ Qi on (t) / R(t) 3 . The damping is caused 
by both, internal heating of the quasi-free electrons in 
the cloud and energy transfer to bound electrons. These 
two effects are responsible for outer and inner ionization, 
respectively. 

For periodic driving F{t) — Fq cos(wi)_the dynamics is 
given by X(t) = A t cos(w£ — (f>t) with 



4_t: 

.18] 



At = -^ 2 ) 2 +( 2r » 2 > 

(fit = arctan(2r^/(f2 t 2 - uj 2 )) . 



(4a) 
(4b) 



The energy balance of the dynamics (@J is characterized, 
on one hand, by energy loss E\ QSS due to the damping 
and, on the other hand, by energy gain -E ga i n from the 
external laser field. The cycle-averaged energy transfer 
rates read 



(E) = (E loss ) + (£ gain ) 
1 
2 



= — Tt At 2 oj 2 + —Fq A t uj sin ( 



(5) 



Obviously and well known |T3 |. maximum (E gain ) or 
optimal heating requires (fit — 7r/2, i. e. resonant be- 
haviour fl t = ui. As shown in Fig. [2Ji, the phase shift 
(fit changes in time from to tt, thereby passing the res- 
onance <fi t = tt/2. This change is directly connected with 
the increased ionization of the cluster, cf. Fig.[2^. In par- 
ticular, the resonance time t rcs coincides with the time 
<2 of maximal ionization rate, cf. Eq. J2J. This applies 
to other laser wavelengths as well (For A — 1170nm t 2 
is somewhat smaller due to the early laser switch-off at 
t = 80 fs.): 
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FIG. 3: Driven and damped harmonic oscillator: amplitude A 
and phase shift as a function of the ratio of eigenfrequency 
Q and driving frequency u> according to Eqs. Q for different 
damping strengths F. 



Passing through the resonance is not necessarily con- 
nected with large amplitude oscillations, if the damping 
strength T t is comparable to or larger than the eigenfre- 
quency fit, see Fig. |3J Otherwise, the phase shift <j>t at 
resonance is independent of the damping (Fig.UJ- 

In order to validate the applicability of the driven and 
damped harmonic oscillator model of Eq. @ we use am- 
plitude A t and phase shift (fit of the collective oscillation 
from our results of Fig. [2] to determine frequency and 
damping rate according to 



£l t 2 = uj 2 + {Fq /A t ) cos (fit, 
T t = (F /(2A t uj))sm<fi t . 



(6a) 
(6b) 



Figure 01 shows these parameters as a function of time 
t for the same system as in Fig. [21 The calculated 
eigenfrequency £lt (gray circles in Fig. 0J) closely resem- 
bles the frequency of a uniformly charged sphere (solid 
line). Moreover, both frequencies match the laser fre- 
quency ui at the same time t f=a 20 fs providing addi- 
tional support for the collective oscillator model. Along 
with the decrease of the eigenfrequency £lt the damp- 
ing term T t (white diamonds in Fig. ^ rises for times 
t w — 60...30fs. This accounts for energy transfer to 
deeper and deeper bound electrons, which does not oc- 
cur in other theoretical studies where also resonant be- 
haviour was discussed 0, [2(j- Either inner ionization 
was not considered or deuterium clusters composed 
of single electron atoms were discussed [2(j. However, 
in order to understand the experimentally observed high 
charge states 0, 13- HI El it is of utmost importance to 
take this into consideration. It is just this continuous 
cycle of effective heating and induced inner ionization at 
the resonance which causes the increased ionization rate 
and the high charge states of the fragments. At this point 
it is interesting to note, that the damping almost com- 
pletely compensates the heating as can be deduced from 
the almost negligible increase of the amplitude of v CM 
before reaching the resonance (Fig. |2t). Once the clus- 
ter has passed the resonance, however, the damping falls 
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FIG. 4: Parameters of the harmonic oscillator model as 
calculated from the Xe923 dynamics in Fig. Solid line: 
eigenfrequency for a spherical, uniformly charged cluster Qt = 
\/Qion(t)/ R(t) a . Circles: eigenfrequency fit l|6a[l . Diamonds: 
damping rate T t d6bl) . Dotted line: laser frequency 10. 

off rapidly (t > 30 fs in Fig. BJ. Obviously, the oscillat- 
ing electron cloud becomes unable to drive further inner 
ionization. This implies a weaker damping and thus an 
increase of the oscillation amplitude, cf. Fig. 

In summary, we have shown that the electron emission 
in medium-sized rare gas clusters (~ 10 3 atoms) is en- 
hanced by resonant energy absorption in agreement with 
experimental data 0, (The metallic nature of the 
clusters used in [i| should be of minor importance for 



the creation of the high charge states > 8 since the delo- 
calized valence electrons are emitted early in the pulse.) 
Our microscopic calculations of the motion of ions and 
electrons using a hierarchical tree code reveal a laser- 
driven collective oscillation of the cloud of quasi- free elec- 
trons which are held back inside the cluster volume by 
the space charge of the cluster. The eigenfrequency of 
this oscillation is determined by charge and size of the 
cluster. Electron emission and cluster expansion change 
these quantities in time and eventually enable the match- 
ing of eigenfrequency Q t and laser frequency uj during the 
pulse. This resonance allows for an effective energy trans- 
fer to the collective motion. The changing phase shift be- 
tween driving field and driven electron cloud clearly indi- 
cates the different stages of energy absorption of the clus- 
ter electrons from the laser field. The fact that the collec- 
tive electron dynamics can be well described by a simple 
damped harmonic oscillator helps to clarify the nature 
of this dynamics and provides a clear signature of this 
type of collective dynamics. We expect that the other 
two mechanisms for effective energy absorption, namely 
enhanced ionization and nanoplasma excitation will have 
an almost vanishing amplitude for the electronic CM ve- 
locity. Small clusters, which exhibit enhanced ionization, 
produce only a few quasi-free electrons which cannot cre- 
ate a sizable CM-velocity amplitude. In large clusters, 
however, where a nanoplasma is formed, electrons are 
heated resonantly at their critical density leaving their 
CM position at rest. 
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